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Dijet production, collision centrality and backgrounds in high-energy p-p collisions 
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Background: Two aspects of high-energy p-p collisions share common phenomenological ele- 
ments: (a) A correlation between jet production and p-p centrality is suggested by the transverse 
partonic structure of hadrons inferred from deep-inelastic scattering data, (b) The underlying event 
(UE) is defined as the final-state particles complementary to a triggered high-energy dijet. An ob- 
servable common to both topics is variation of so-called transverse multiplicity N± with SiPt,trig dijet 
trigger. Purpose: Test assumptions associated with p-p collision centrality and the UE. Determine 
the nature of the UE and explore the relation between jet production and p-p centrality. Method: 
Use the two-component model (TCM) of spectra and correlations derived from 200 GeV p-p collisions 
to construct a simulated particle distribution on {pt,nch)- Use the simulation to predict the N±_ 
response to pt,trig- Use the measured pt spectrum Uch dependence to explore the effect of changing 
p-p centrality. Results: The TCM provides a good description of N± vs pt,trig and the N±{pt) 
spectrum. The relation of N± to pt and pt,trig departs from assumptions. Conclusions: The pt 
spectrum TCM combined in this analysis with measured minimum-bias p-p angular correlations 
suggests that the UE includes a substantial contribution from the triggered dijet in addition to 
the contribution from projectile fragmentation (beam-beam remnants). The jet contribution to N± 
may represent a universal large-angle base common to all dijets that extends across 2n azimuth. 
The analysis further suggests that p-p centrality is not controlled significantly by pt,trig but may 
be correlated to some extent with an imposed n^h condition, depending on the role of fiuctuations. 
Future correlation studies may better determine the role of p-p centrality. These results may have 
implications for ongoing RHIC analysis and LHC searches for physics beyond the standard model. 

PACS numbers: 12.38.Qk, 13.87.Fh, 25.75.Ag, 25.75.Bh, 25.75.Ld, 25.75.Nq 



I. INTRODUCTION 

The hadronic final state of high-energy p-p collisions is 
of fundamental importance as a reference for interpreta- 
tion of nucleus- nucleus collisions at the RHIC and LHC. 
The search for novel physics in any new collision system 
or energy regime must refer to what is "conventional" in a 
standard-model QCD description of elementary hadronic 
collisions. There are several possible contributions to 
single-particle spectra and two-particle correlations, in- 
cluding transverse fragmentation of large-angle-scattered 
partons to jets, longitudinal projectile-nucleon dissocia- 
tion, radial and other flow components, Bose-Einstein 
correlations (BEC) and possibly other processes not yet 
anticipated theoretically. 

Of particular interest are two possibly-related issues: 
(a) Are there significant manifestations of p-p collision 
centrality in final-state hadron spectra and correlations? 
Is the concept of collision centrality viable in a small 
system possibly dominated by fluctuations? If so, how is 
p-p centrality related to final-state observables, especially 
jet production? (b) What is the nature of the underlying 
event (UE) , those p-p collision products characterized by 
transverse multiplicity N±_ (relative to the dijet axis) that 
are assumed to be distinct from a triggered high-energy 
dijet? What physical mechanisms may contribute? 

In Rcf. Pi deep-inelastic scattering (DIS) data are in- 
terpreted to suggest that the low-x giuon component of 
the proton wave function is confined to a radius substan- 
tially smaller than the overall hadron size. Jet production 
at mid rapidity in p-p collisions might then correspond to 



reduced impact parameter. A pt trigger condition pt,trig 
should favor jet production at higher pt- A nonjet com- 
ponent of the event multiplicity, transverse multiplicity 
N± , could be a measure of collision centrality. A specific 
trend of N± vs pt.trig might then confirm the conjectured 
centrality dependence and identify the pt,trig interval fa- 
voring jet production. 

A high-pt trigger condition may serve as a form of jet 
selection where event-wise jet reconstruction is imprac- 
tical. But other methods are available to isolate or in- 
fer jet systcmatics, including differential spectrum analy- 
sis |2jj3j and combinatoric two-particle correlation anal- 
ysis [jj^l ■ In Ref. a two-component (soft -I- hard) anal- 
ysis of Pt spectra for several p-p event-multiplicity classes 
revealed a spectrum "hard" component that scales with 
multiplicity as n^,j and has properties consistent with a 
minimum-bias jet contribution (no condition on jet en- 
ergy) . The "soft" component scales linearly with rich and 
is consistent with projectile-proton dissociation (beam- 
beam remnants). The relation of the hard component to 
perturbative QCD (pQCD) predictions was established 
in Ref. Two-particle correlations from p-p collisions 
have also been related to expected jet systematics 4, §]. 

In the present study we consider the proposal from 
Ref. [1] in the context of the pt spectrum systematics de- 
scribed in Ref. ^] . We use spectrum information to gen- 
erate jet-related (hard) and nonjet (soft) contributions to 
hadron production distributed on the 2D space (j>t,nch)- 
We then determine the effect of a pt,trig condition on 
soft and hard components compared to expectations from 
Ref. [l[. The present analysis could be extended to jet 
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structure in 2D angular correlations [1, |^ 0| and a non- 
jet quadrupole component observed to depend strongly 
on collision geometry in A- A collisions [1, Q . 

This study also addresses the underlying event. Varia- 
tion of UE properties (azimuth correlation structure and 
A^j^) with event selection based on a pt,trig condition has 
been studied extensively |lCll- ll2| . There is a formal rela- 
tion with the proposal of Ref. [l| in terms of N± response 
to variation of pt,trig- In both cases N± is assumed to 
be unrelated to the triggered dijet and sensitive to the 
properties of the underlying event. In the present anal- 
ysis we combine the event-type selection model devel- 
oped in relation to Ref. [l| with two-component spec- 
trum and angular-correlation phenomenology from 200 
GeV p-p collisions [3, ^ to develop an alternative inter- 
pretation of N± systematics and the nature of the UE. 

Finally, the measured Uch systematics of jet production 
in spectra and correlations from 200 GeV p-p collisions 
are compared with a pQCD description of dijet produc- 
tion in the context of a Glauber model based on trans- 
verse parton distributions within projectile nucleons. We 
seek evidence from the trend of jet production vs Uch for 
some relation to p-p centrality. The event-wise condi- 
tion is complementary to the pt,trig condition employed 
in UE studies and proposed in Ref. [l| to demonstrate 
the relevance of p-p centrality to jet production. 

This paper is arranged as follows: In Sec. |lT]we intro- 
duce two descriptions of jet production in p-p collisions 
in the context of proton transverse structure, p-p central- 
ity and the underlying event that provide a context for 
the present study. In Sec. lIIII we review a two-component 
model of particle yields, spectra and correlations inferred 
from 200 GeV p-p collision data. In Sec. IIVI we ana- 
lyze the structure of the transverse multiplicity spectrum 
dN±{pt)/dpt for triggered events at 0.9 and 1.8 TeV. In 
Sec. |V] we apply the two-component model to the p-p 
minimum-bias event distribution on Uch- 

In Sec. IVII we construct ensemble-averaged particle 
densities on event multiplicity Uch and transverse rapid- 
ity Tjt for soft and hard hadron production mechanisms 
and for soft (no jets) and hard (at least one jet) event 
types. In Sec. IVIII we present measured angular corre- 
lations for minimum-bias jets in 200 GeV p-p collisions, 
emphasizing their contribution to the transverse azimuth 
region and N±. In Sec. IVIIII we illustrate the relation 
between transverse multiplicity A^_l and trigger pt cut 
that is a central element of underlying-event studies. In 
Sec. IIXI we review the Glauber model applied to Au-Au 
and p-p collisions and its implications for a correlation 
between particle and dijet production and p-p collision 
centrality. In Sec.|X]we compare the results of this study 
to the material presented in SecHT] Sec. IXII summarizes. 



II. BASIC STRUCTURE OF p-p COLLISIONS 

We consider descriptions of p-p dynamics and final- 
state structure in two contexts. A theoretical description 



is based on deep-inelastic scattering (DIS) data and con- 
siders possible manifestations of p-p collision centrality in 
certain observables. An experimental description is based 
on the azimuth correlation structure of events containing 
a triggered or reconstructed dijet, especially properties 
of the underlying event assumed to be complementary to 
the dijet. Each approach relies on a pt.trig condition. 

A. A theoretical context 

In Ref. [H it is argued that hard particle production 
(i.e., large-angle gluon scattering to dijets) should be 
most probable at small p-p impact parameter because 
the transverse size of the low- a; gluon distribution in the 
proton inferred from DIS data is substantially smaller 
than the overall proton size. Soft particle production 
(not associated with a triggered dijet) should vary with 
b over a large range and may serve as a centrality mea- 
sure. Specifically, transverse multiplicity N± (perpendic- 
ular to a trigger-particle momentum and therefore to a 
dijet axis) may be correlated with b. Absent a jet trigger 
collisions should be dominated by soft processes at larger 
6 with N± Nnsd (NSD denotes non-single-diffractive 
p-p collisions). 

It is assumed that jet production should be correlated 
with smaller b and therefore larger N±. The larger mul- 
tiplicity may then make it difficult to observe lowev-pt 
jets directly, and semihard dynamics, and particle pro- 
duction with Pt « few GeV is said to be not well un- 
derstood. Therefore, indirect selection of hard processes 
(jets) is established with a pt condition denoted here by 

Pt,trig- 

For sufficiently high pt.trig b should be relatively small 
and nearly independent of the trigger condition. Equiva- 
lently, transverse multiplicity N± should be nearly inde- 
pendent of pt,trig and substantially larger than for NSD 
p-p collisions dominated by soft proton dissociation at 
large b. Reference [l| then poses the question (given the 
several assumptions): above what critical pt,trig value is 
hadron production dominated by "hard" parton-parton 
interactions in more-central p-p collisions? In the present 
analysis we examine those assumptions and expectations 
in the context of differential analysis of pt spectra vs p-p 
multiplicity and the correspondence with established jet 
phenomenology and pQCD predictions. 



B. An experimental context 

The CDF collaboration has emphasized study of the 
underlying event in p-p collisions [l^, [IH. In colli- 
sions that include a (triggered) hard parton scatter to 
a high-energy dijet the UE is defined as what is not 
part of the dijet, with several possible production mech- 
anisms including the beam-beam remnant (soft compo- 
nent, projectile proton dissociation), multiple parton in- 
teractions (MPI, untriggered semihard parton scatters. 
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see Ref. [13| for a review) and initial-state radiation. The 
UE is explored in the transverse region (TR) said to be 
"very sensitive" to the UE [l^. The TR, defined by 
— 7r/2| < tt/6, includes total multiplicity N± within 
some 77 acceptance A77 and is assumed to be outside the 
triggered hard-scatter dijet on azimuth. Two other com- 
plementary azimuth regions are denoted toward (includ- 
ing the same-side or triggered jet) and away (including 
the away-side or recoil jet). 

CDF UE analysis is based on a particle acceptance 
defined by pt > 0.5 GeV/c and I77I < 1. Jets are defined 
by a cone jet finder with radius i? < 1 on (rj, (p) and 
include jet fragment multiplicities as low as 1. Relevant 
issues for UE analysis include particle production and 
correlation mechanisms, trends of N±_ on collision energy 
and trigger pt condition and the iVj^ pt spectrum. 

Some results from the CDF UE analysis at ^/s = 1.8 
TeV are particularly relevant to the present study; In 
Fig. 13 of Ref. the growth of total event multiplicity 
within the ry acceptance is plotted vs pt,trig- In Fig. 21 
the total yield is separated into three parts on azimuth 
relative to the trigger. The N± yield from the transverse 
region assumed to represent the UE increases with pt,trig 
from zero to a limiting plateau value ~ 2.3. The increase 
of N±{pt,trig) is attributed to an MPI rate increasing 
to a saturation value. In Fig. 17 the azimuth distribu- 
tion is shown for a specific jet trigger [Pj > 5 GeV/c) 
in relation to the three equal azimuth regions. And in 
Fig. 37 the N± pt spectrum for the TR is shown for three 
trigger conditions. The N± spectrum is assumed to be 
monolithic (described by a single functional form) . In the 
present study we consider those results in the context of 
minimum-bias systematics from 200 GeV p-p collisions 
and a two-component model of hadron production. 



C. Strategy for the present analysis 

We first present two-component models for p-p pt 
spectra {rich dependence), angular correlations (pt de- 
pendence) and minimum-bias event distributions on rich 
[based on inferred event- wise jet number rij^rich)]- From 
those elements we construct particle-number densities on 
the space (pt,rich) (particle transverse momentum and 
event multiplicity) that form the basis for studying rela- 
tions among a,pt,trig condition, particle production mech- 
anisms and event types. For a given pt^trig value what 
are the corresponding fractions of particles from soft and 
hard production mechanisms, and the fractions of soft 
and hard event types? How are those trends related 
to the systematics of transverse multiplicity N±7 We 
also consider the consequences of a complementary event- 
multiplicity rich condition on spectra and correlations in 
the context of a Glauber model for pQCD applied to p-p 
collisions. 



III. THE TWO-COMPONENT MODEL 

The two-component model of p-p pt spectra was ini- 
tially derived from the multiplicity dependence of spec- 
tra [31 • Its physical interpretation is based on a sim- 
ple model of hadron production: longitudinal projectile- 
nucleon dissociation (soft) and transverse large-angle- 
scattered parton fragmentation (hard). In more- 
peripheral A-A collisions hadron production from the two 
processes scales as A^part/2 (number of nucleon partici- 
pant pairs) and N^in (number of binary N-N collisions) 
respectively 0, • Those parameters are related to ob- 
servable rich via the Glauber model of A-A geometry jlSi] . 
For p-p collisions the corresponding geometry parame- 
ters and their relation to observables are considered in 
Sec. IIXI Interpretation of model elements in terms of 
physical mechanisms is fairly reliable for p-p collisions 
but is less certain for more-central Au-Au collisions. The 
TCM can also be applied to single-particle yields and 
two-particle correlations. We require that all TCM in- 
terpretations be consistent across spectra, correlations, 
minimum-bias distributions on rich and particle yields. 



A. p-p Pt spectrum systematics vs nch 

In Ref. ^ the structure of pt spectra for several charge- 
multiplicity Uch classes of p-p collision events was ana- 
lyzed. The spectrum data are well-described by a two- 
component model function of the form 

duch/ytdyt = ris{nch)So{yt) +nh{nch)HQ{yt) (1) 

(within rapidity acceptance Ar; = 1), where unit-integral 
factors So{yt) and Ho{yt) (defined in App.|X| are respec- 
tively a Levy distribution on mt (suitably transformed 
to yt) and a peaked distribution on yt approximated 
for this analysis by a Gaussian,. Transverse rapidity 
yt — ln{(mt -\-pt)/mh} is adopted to aid visualization 
of low-pt spectrum structure, and pion mass ruh ^ m^r 
is assumed for unidentified hadrons. Soft rig and hard 
Uh multiplicities observed within Arj = 1 approximate 
densities drix/drj. Quantities Hpp{yt) = nhHolyt) and 
Spp{yt) = risSoiyt) are also defined. Equation ([T]) repre- 
sents approximate factorization of the multiplicity and yt 
dependence, effectively in the form of a two-term Taylor 
expansion on rich- Series coefficients for higher terms as 
determined by spectrum data are consistent with zero. 

Figure [T] (left panel) shows the most important re- 
sult from Ref. [2]. Data in the left panel were obtained 
by dividing each pt spectrum by the soft-component 
yield risirich) (inferred iteratively) and subtracting the 
fixed soft-component function So{yt) derived in the limit 
"n-ch ^ 0. The result is inferred hard components Hpp/ng. 
The data are consistent with rih/ris oc rich and rih n^. 

Figure [1] (right panel) shows a combination of data 
from the left panel corresponding to non-single-diffractive 
(NSD) p-p collisions (points). The solid curve through 
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FIG. 1: Left: Evolution of the pt spectrum hard component 
with event multiplicity rich, from 200 GeV p-p collisions 
Right: Calculated pQCD fragment distribution (solid curve) 
corresponding to the NSD-equivalent hard component from 
the left panel (points) 0]. The dash-dotted curve is the cor- 
responding Gaussian form from the left panel (solid curves). 



the data is derived by combining published fragmentation 
functions |l3 with a pQCD parton energy spectrum @. 
The dash-dotted curve is the simplified Gaussian form 
from the left panel, defined in App. |X] and used in this 
analysis. The good agreement between pQCD theory 
and p-p spectrum data strongly supports interpretation 
of the spectrum hard component as a jet fragment distri- 
bution representing a minimum-bias parton energy spec- 
trum bounded below at 3 GeV. We conclude from Ref. [H 
that varying the p-p Uch condition controls the jet fre- 
quency. Jet number nj is proportional to n^nch ~ '^l.h^ 
and the constant of proportionality is consistent with 
pQGD and angular acceptance Ary = 1. 



B. p-p two-particle correlations 



Extensive measurements of p-p two-particle angular 
correlations on (ry, 0) and correlations on transverse ra- 
pidity space (?/ti,2/t2_) for several charge combinations 
have been reported 0, Q • Most of the resulting struc- 
ture is well-described by the TCM. The exceptions are 
contributions from Bose-Einstein correlations and elec- 
tron pairs from photon conversions. The soft component 
appears mainly below pt = 0.5 GeV/c whereas the hard 
component appears only above 0.35 GeV/c. Thus, a cut 
at 0.5 GeV/c imposed on both particles of each pair ef- 
fectively separates hard and soft components. The hard 
component on (1/41,2/42) projected onto ID yt is compat- 
ible with the structure in Fig. [TJ and the corresponding 
shapes on (77, (/>) are as expected for interjet and intrajet 
angular correlations, strongly supporting a jet interpre- 
tation for the correlation hard component, p-p angular 
correlation data are discussed further in Sec. I VIII 



C. Two-component model and A-A collisions 

The present study emphasizes jet production and the 
role of collision centrality in p-p collisions. Centrality 
trends in A-A collisions have been studied extensively. 
The relevant centrality observable for A-A collisions is 
Uch, and the corresponding Glauber centrality measure 
is mean participant path length v = 2Nbin/Npart- The 
two parameters are related by the total cross section [l^ . 
By analogy with A-A collisions Uch may also be the pre- 
ferred control parameter in p-p collisions if centrality is a 
relevant concept there. Whereas dijet production is ob- 
served to scale as n^^ in p-p collisions it scales as Nmu 
in more-peripheral A-A collisions (as expected for lin- 
ear superposition of N-N collisions). From the Glauber 

model we obtain Nhin ~ A"pirt ^nd Npart Uch (ap- 
proximately). Therefore dijet production scales as ri^j^ 
in more-peripheral A-A collisions, not n^f^ as in p-p col- 
lisions. Scaling issues are discussed further in Sec. IIXDI 



D. Two-component theory description 

Reference [l| introduces an alternative two-component 
model for transverse multiplicity N±_ systematics based 
on a two-state system of large-6 events with no jets (soft) 
and small-6 events with jets (hard). The form is 

N±{pt,trig) = K{Pt.,trig)Ns{hs) + )^h{Pt,trig)Nh{bh){'^) 

with fractions A/j + Ag = 1 of hard and soft events respec- 
tively for a given pt,trig and fixed numbers Ns and Nh, 
with Nh ^ Ns assumed. While the form of that equa- 
tion is similar to Eq. ^ the underlying assumptions are 
quite different. By definition of N±{pt^trig) (complemen- 
tary to the triggered dijet) both Ns (larger b) and N^ 
(smaller b) include soft (nonjet) components from the 
two event types. An MPI contribution to Nh is also pos- 
sible. The subscripts refer to the event types with their 
different centralities, not soft and hard components of a 
given event. Fractional abundances A3; of the two event 
types are assumed to be determined only by pt,trig- 

In the present analysis we define two types of p-p col- 
lision: soft and hard events. Equation ([!]) applies to any 
p-p event. Hard events contain at least one minimum-bias 
jet within some specified detector angular acceptance. 
Soft events contain no jet structure within the accep- 
tance {ufi 0). All events include a common nonjet soft 
component of particle production. The relation of soft 
and hard events to p-p collision centrality (if any) is not 
assumed a priori. We compare theory assumptions with 
data trends and find that the soft components of N± for 
the two event types are actually similar, but N± for hard 
events includes an additional triggered-jet-Telated (hard) 
component that may even dominate N±. 
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IV. TRANSVERSE MULTIPLICITY Nx(pt) 

We now apply the two-component model to study 
transverse multiplicity N±(pt) spectrum structure and 
other systematic issues. We assume that measured pro- 
jectile dissociation (beam-beam remnants) contributes to 
the N± spectrum according to the TCM but also con- 
sider a possible hard-component contribution that may 
arise from the triggered dijet. We find that the TCM 
provides an accurate description of N±{pt) spectra and 
leads to an alternative interpretation of UE composition. 

A. Nx(pt) spectrum structure 

Figure [5] shows transverse multiplicity N±{pt) data for 
1.8 TeV (solid points) for a specific trigger condition 
(summed jet momentum Pj > 5 GeV/c) from Fig. 37 
of Ref. [13 ■ The data can be described with TCM 
Eq. ([T]) modulo transformation yt — >■ pt and additional 
factor Pt- The curves are unit-normal forms pt So{pt) 
and Pt Ho{pt) derived in Ref. from -y/s = 200 GeV p-p 
data. The respective coefficients are N±,s — 3.3 ~ 2x5/3 
(dashed curve) and N±^h = 1-3 ~ 4/3 (solid curve). The 
sum (dash-dotted curve) describes the 1.8 TeV data well. 
Data for 0.9 TeV with the lower 3 GeV/c trigger condi- 
tion (open points) are discussed in Sec. IX Al 

The Eq. ^ coefficients (with systematic uncertainties) 
were determined from the data as follows. The spectrum 
shapes S'nfpt) and -ffo(Pt) inferred from the 200 GeV data 
m Ref. (see App. El were left unchanged. Coeffi- 
cient N±^h 1-3 was adjusted to fit the CDF data (sofid 
points) at larger pt, implying that n;i = 3x 1.3 « 4±0.5 is 
the hard-component fragment integral within 27r azimuth 
and A?7 = 2, interpreted as the minimum jet fragment 
multiplicity for imposed conditions (e.g., at least one jet). 
Coefficient N±,s — 3.3 « 2 x 5/3 was then adjusted to 
fit the data at smaller pt dominated by the lowest point. 
The corresponding nonjet soft density dus/drj « 5 ± 1 
can be compared with the NSD soft multiplicity density 
dus/dr] w 4 measured at 1.8 TeV [lo| . 

Figure [5] (right panel) shows the same data and curves 
on transverse rapidity yt with factor 1/yt added. While 
Pt may be a directly measured quantity it obscures details 
at smaller pt as compared with yt . The linear power-law 
trend evident at larger yt (compare with the dotted line 
at right) is consistent with an underlying minimum-bias 
pQCD parton spectrum @ . 

The dash-dotted curve in the right panel can be com- 
pared with Fig. 10 (right panel) of Ref. predicting the 
TCM spectrum shape for 200 GeV p-p hard events, those 
with at least one jet (dotted curve). The soft component 
at 1.8 TeV is increased relative to 200 GeV by 5/2.5 = 2 
times whereas the minimum-bias hard component (jets) 
may be more slowly varying with collision energy. Al- 
though the minimum in the hard-event spectrum from 
Ref. 01 does not appear in Fig. [5] (right panel) the two 
results are mutually consistent within uncertainties. 
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FIG. 2: Left: Transverse multiplicity N_i_ spectrum on pt for 
1.8 TeV p-p collisions (solid points, I77I < 1) from Fig. 37 of 
Ref. [lO] and for 0.9 TeV p-p collisions (open points, \rj\ < 2) 
from Fig. 4 of Ref. [l2|. Curve 5* is the soft component with 
data normalization corresponding to a soft-component density 
dus/dri = 5. Curve H is derived from the pQCD-calculated 
fragment distribution for minimum-bias jets from 200 GeV 
p-p collisions 0| also normalized to match the 1.8 TeV data. 
Right: Data and curves from the left panel transformed to 
transverse rapidity yt with the appropriate Jacobian. A factor 
1/yt is added to compare with Fig. [T] (right) and Eq. {T]). 



B. Collision energy comparisons 

In Fig. [5] (right panel) the vertical scale range is shifted 
by a factor 30 higher than that for Fig. [T] (right panel) 
to explore the relation between the two data sets. The 
jet frequency inferred for the untriggered 200 GeV data 
in Fig. [H is / = 0.013 [1], and the spectrum is inte- 
grated over 27r azimuth and A77 — 1. The estimated 
most-probable jet energy is 3 GeV. The triggered 1.8 
TeV data in Fig. [5] with trigger-jet momentum Pj > 5 
GeV/c correspond to jet frequency / ~ 1 and are inte- 
grated over 27r/3 azimuth and A77 = 2. Setting aside the 
difference in A77 the nominal relation between data sets 
is then 1/(3 x 0.013) « 25 The functional form de- 
rived from pQCD (sohd curve, ^) appears to describe 
the TR hard-component spectrum accurately. 

The hard-component maximum value in Fig. [1] (right 
panel) is about 0.01 whereas in Fig. [5] (right panel) it is 
about 0.4, with identical shapes. The ratio 40 should be 
compared with the factor 25 derived above. The hard- 
component yield at 200 GeV is apparently about 60% 
of that at 1.8 TeV. But there are two related issues: (a) 
effects of different 77 acceptances and (b) bias effects from 
a 5 GeV/c jet trigger vs no trigger. The apparent factor- 
1.6 difference may simply arise from differences in the jet 
rj acceptance and effective mean jet energy. Those issues 
are discussed further in Sees. IVII^ and IVini 



C. Spectrum extrapolation 

N±{pt) spectrum extrapolation in Sec. V of Ref. [l3] 
assumes a monolithic spectrum approximated by 
exp(— 2pt) (dotted curves in Fig. [5]). For plateau value 
N± = 2.3 the inferred 77 density is dnch/di] = 3.5 with 
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Pt > 0.5 GeV/c. Extrapolation to pt = with the as- 
sumed spectrum shape gives drich/dri « 10 compared 
with soft component dus/dri w 4. The ratio 2.5 suggests 
a substantial unexplained contribution to the UE ilO]. 

However, accurate extrapolation of any pt spectrum 
requires knowledge of its structure. Measured p-p spec- 
tra at 200 GeV can be decomposed into soft and hard 
components with known spectrum shapes, as in App. |^ 
We apply the same TCM to N±{pt) assuming the com- 
ponent shapes are similar at 1.8 TeV. As noted, the soft- 
component coefficient derived from the spectrum data in 
Fig. [5] corresponds to dng/dij « 10/2 = 5. The hard- 
component coefficient N±ji corresponds to duh/drj sa 
4/2 = 2. The integral of the apparently unexplained 
(hard) part of the UE is thus only 2/7 « 30% of the to- 
tal, and according to the TCM that component may be 
identified in some sense with jet production. 

The most surprising aspect of the UE pt spectrum is 
not a factor-2.5 ratio of extrapolated N±_ to that expected 
for soft events. It is the factor-20 ratio at yt ~ 3.3 [pt ~ 2 
GeV/c) between the N±{pt) spectrum (solid curves) and 
the expected soft spectrum component (dashed curves) in 
Fig. [51 That result strongly suggests that a large jet- 
related (hard) contribution may dominate the UE. The 
origin of the hard contribution to N± (MPI or triggered 
dijet) is a central issue for this study. 

V. p-p MINIMUM-BIAS DISTRIBUTION 

Equation ([1]) describes the pt spectrum structure for 
a given event multiplicity rich- We also require a model 
that describes the minimum-bias event distribution on 
rich for 200 GeV p-p collisions and its soft and hard com- 
ponents. We adopt a model consistent with measured 
energy systematics up to 7 TeV. We then use the ob- 
served jet production rate vs rich from Ref. 0] to deter- 
mine the fraction of soft events Po{nch) (no jet in the 
acceptance) and of hard events 1 — Po{nch) (at least one 
jet in the acceptance), with Po{rich) — exp[— nj(ncft)] and 
nj{nch) the mean jet number in angular acceptance A77 
for event multiplicity rich as reported in Sec. lIX Cl Quan- 
tity Nt = cTxCt represents the total number of events in 
an ensemble for some trigger condition X (e.g., NSD) and 
integrated luminosity Ct. The unit-normal event distri- 
bution on rich is represented by {\/Nt)dNt/dnch- 

Figure [3] (left panel) shows distributions on rich/ ^r] 
from the CMS collaboration (broken curves) for three 
collision energies [3| • The data approach asymptotically 
the solid curve estimating the minimum-bias distribution 
for 200 GeV collisions. The 200 GeV parametrization is 

^^ = 0.16exp(-|n,„/A,7-2.3|iV5). (3) 

JVt drich 

Although an approximation to the data, the model func- 
tion is convenient and suffices for the present study. Fig- 
ure [3] (right panel) shows an expanded view of the 200 
GeV model with the soft and hard components. Ns is the 




FIG. 3: Left: Measured minimum-bias distributions on rich 
for three energies (broken curves) from Ref. with \r]\ < 
2.4 and an extrapolation to 200 GeV (solid curve). Right: 
The 200 GeV extrapolation T is separated into soft S and 
hard H components based on results from Ref. [2[ (event- 
wise jet frequency /) providing a two-component model of 
the minimum-bias distribution. 

number of soft events (distribution S), and Nh the num- 
ber of hard events (distribution H) with Nt — Ng + Nh 
based on Poirich) defined above. 

The format of Fig. |3] with accepted rich normalized by 
A77 contrasts with that in Ref. jl8| where total multi- 
plicities from different 77 acceptance intervals are com- 
pared. Use of rich I A77 (consistently on both abscissa and 
ordinate) may aid comparison of data from different ac- 
ceptances and energies. A range of 77 acceptances (e.g., 
A77 = 1 4.8 in Ref. [l^) implies nearly a factor- two 
variation in the number of jets per detected dijet (i.e., the 
probability that a jet within the 77 acceptance also has the 
recoil partner within the acceptance) [l^l . Figure [3] can 
be compared with the so-called KNO format [1^ wherein 
nPn is plotted vs njfi [P„ is here [1/ Nt)dNt/ drich]- The 
acceptance A?7 does cancel in both product and ratio, 
but the mean rich may include contributions from two or 
more production mechanisms, making overall interpreta- 
tion difficult. 



VI. TOTAL PARTICLE DENSITY ON (yt.nch) 

Wc have defined two-component models for particle 
multiplicity distributed on yt as spectra drich /ytdyt, and 
for event number distributed on event multiplicity as 
dNtI drich- We now combine those results to produce a 
2D distribution of total particle multiplicity Nch = fichNt 
(defining ensemble- mean event multiplicity rich) on the 
space {y It rich)- The 2D particle distribution can be de- 
composed into soft and hard components in two ways: 
according to hadron production mechanism (particle ori- 
gin) and according to event type. A distinction should 
be maintained between soft and hard particle production 
processes and soft and hard event types. A hadron from 
a soft process may emerge from a hard event, but the 
reverse is not possible. 

In the following we assume that a trigger particle for 
given yt,trig bin is represented on average by the most- 
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probable particle in that yt bin from all events in the en- 
semble, which then defines the most-probable event rich 
(the mode) for that condition. Those events are most 
likely to produce the required trigger particle. 



A. Description according to particle type 

Decomposition according to particle type (hadron pro- 
duction mechanism) can be expressed by 



dNf dn. 



ch 



duch Vtdyt 



dNt 

duch 
dNt 

dn, 



ch 



ns{nch)So{yt) 
nh{nch)HQ{yt), 



(4) 



where a common factor Arj/Nt is suppressed, — 
nch — nu and Uh — nj{nch) nch,j is the hard component 
of event multiplicity (jet fragments), drich appearing in 
a denominator and referring to an event index should be 
distinguished from drich in a numerator denoting a dif- 
ferential element of event multiplicity. 

Figure m shows surface plots of the soft component (left 
panel) and hard component (right panel) of the total mul- 
tiplicity distribution on {yt, rich) defined by Eq. (jj]) (first 
and second terms respectively). The z-axis scales are 
logarithmic. Dashed lines indicate soft and hard distri- 
bution modes (most probable values) nch,s/Av — 3-5 and 
nch,h/Av = 5.7. A condition imposed on yt would con- 
trol the distribution shape on rich, whereas a condition 
on rich would control the yt spectrum shape. We want 
to determine what event class (e.g., rich value) is most 
probable for a specific trigger condition on pt or yt . 
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FIG. 4: (Color online) Left: 2D surface representing the 
Cartesian product of two ID distributions, the first term (soft) 
ofEq. Q. Right; Surface representing the second term (hard) 
of Eq. (Q. The dashed lines represent the modes on rich I ^t] 
at 3.5 and 5.7. 

Figure [5] (left panel) shows the sum of the two terms 
in Eq. The vertical dashed lines indicate the modes 
on rich/ A'n of the soft and hard components. It is no- 
table that the separation is approximately 2.5, the mean 
fragment multiplicity of a minimum-bias jet or minijet. 
The dotted lines indicate three yt conditions employed in 
discussion. The bold solid line indicates the CDF accep- 
tance cutoff at Pt =0.5 GeV/c. 
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FIG. 5: (Color online) Left: 2D surface representing the sum 
of soft and hard terms in Eq. Q. The dashed lines indi- 
cate the soft and hard modes on rich /A77. The dotted lines 
represent three yt conditions for illustration. Right: Corre- 
sponding particle densities on event multiplicity nch/^rj for 
the three conditions (cuts) on yt (with pt « 0.9, 1.4 and 2.2 
GeV/c respectively). Solid curves (labeled by yt cut) repre- 
sent total densities. Dashed curves represent soft components 
and dash-dotted represent hard components. 



Figure [5] (right panel) shows conditional particle den- 
sity distributions on event rich as defined by Eq. ^ for 
three conditions (trigger cuts) on yt. For the lowest cut 
the total density is dominated by the soft component 
(dashed). For the highest cut the total is dominated 
by the hard component (dash-dotted). Variation of a 
Pt,trig condition over a rather small interval can change 
the particle-type mixture from almost all soft-component 
hadrons to almost all hard-component hadrons. How- 
ever, the soft component of the event multiplicity (pos- 
sibly related to p-p centrality) does not change signifi- 
cantly. That panel can be compared with the spectrum 
trends in Fig. 1 of Ref. [2| and in Fig. [T] (left panel) of 
this paper where rich/Ar] is the condition variable and 
the soft-component yield varies over a large range. 



B. Description according to event type 

In Ref. m it is proposed to control the fraction of 
events that include hard processes (hard events) by spec- 
ifying a Pt trigger condition. With increasing pt.trig the 
fraction of hard events should increase. Assuming a 
strong correlation among jet production, p-p centrality 
and soft multiplicity there should also be a substantial 
increase of transverse multiplicity N± (assumed compa- 
rable to the multiplicity soft component of the TCM) 
up to some critical pt.trig value. Beyond that value ac- 
cepted collisions should be mainly hard, nearly central 
and therefore with similar characteristics including iVj^. 

In Eq. ([3]) the 2D particle density is decomposed ac- 
cording to hadron production mechanism. Decomposi- 
tion according to event type can be expressed by 

— Po{nch)^-^ ns{nch)So{yt) 5) 

dnchVtdyt Nt drich 

= Cs{nch)So{yt) 
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for soft events (no jets), with 



rich since Uh = 0, and 
(6) 



— — = [I- Poinch)]^-: X 

dnchUtdyt Nt duch 

= C'^{nch)So{yt) + Ch{nch)Ho{yt) 



for hard events (at least one jet), where n(, — rich — n' 
and n'y^ — n?i/[l — -Po(f^c/i)] > ^ch,j is distinct from rih 
averaged over all soft and hard events combined. 
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FIG. 6: Left: Factors from Eqs. ((5| and ((6| that multiply 
and Ho respectively, ps describes the soft component 
of soft events and describes the soft component of hard 
events. Right: Corresponding event-type fractions derived 
from Eq. ((Tjl corresponding to definitions in Ref. ^J. 

Figure [6] (left panel) shows factors Cs{nch), C'g{nch) 
and Chinch) defined by Eqs. ([5]) and For soft events 
we observe that nch,sl At] = 3.2 and the maximum soft- 
event particle density Cs,max = 0.39. For hard events we 
obtain nch,h/ Arj = 5.7 and maximum hard-event particle 
density Ch,max — 0.035. The hard and soft event-type 
modes nch,x are formally comparable to asymptotic val- 
ues Nsoft and Nhard in Ref. 1], but the actual hadron 
production mechanisms and event-type criteria may not 
be as assumed. 

Figure [S] (right panel) combines properties of the multi- 
plicity densities to obtain estimated event-type fractions 
Aa; vs yt,trig as introduced in Eq. ([2]). Based on the TCM 
the soft-event fraction is 



As(y; 



t,trig ) 



Cg.max S{^{yt^trig ) 



So{yt,trig) + Ch .max 



(7) 



with a corresponding expression for Xh{yt,trig) ii s,S <->■ 
h,H. Those results relate to Eq. (15) of Ref. [ll]. Since 
is at most 10% of Cs it is omitted from Eq. (jT]). The 



event- type fractions cross at pt. 



trig 



1.5 GeV/c, and 



above 2.5 GeV/c hard events comprise more than 90% of 
the total. Although only 3 % of 2 00 GeV NSD p-p colli- 
sions are hard events (Sec. IIXCI) a pt,trig > 2.5 GeV/c 
condition selects almost all such events and almost no soft 
events. The soft component of the hard-event ensemble 
is approximately that of the minimum-bias sample. 



VII. ANGULAR CORRELATION STRUCTURE 

We can use the structure of measured minimum-bias 
angular correlations from 200 GeV p-p collisions to pre- 
dict a possible dijet contribution to the transverse az- 
imuth region at higher collision energies, assuming that 
minimum-bias jets are approximately universal over a 
broad range of p-p collision energies. 



A. Minijets from p-p collisions at 200 GeV 

Figure [7] (left panel) shows a parametrization of mea- 
sured 2D angular correlations representing minimum- 
bias jets (minijets) in 200 GeV NSD p-p collisions from 
Refs. 0, No trigger condition is imposed — the distri- 
bution represents all combinatoric pairs above 0.5 GeV/c 
(after subtraction of a uniform background based on 2D 
model fits) . Figure [7] (right panel) shows the 2D his- 
togram projected onto ID azimuth (solid curve). Note 
that the away-side pair structure (dashed curve) has been 
modified by a triangular 7y-acceptance correction that 
overestimates the accepted AS pairs. The dash-dotted 
sum represents the uncorrected (accepted) pair density. 
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FIG. 7: (Color online) Left: Minimum-bias jet-related 2D 
angular correlations from 200 GeV p-p collisions 0, IBl- A 
same-side 2D peak at the origin is elongated on azimuth. The 
away-side ID peak is broad on azimuth and approximated by 
a dipole cos(/!>a form. Right: Projection of jet-related angular 
correlations onto azimuth. The hatched regions represent the 
periodic "tranverse region" (TR) invoked in underlying-event 
studies. The unhatched regions, denoted "toward" ((^a ~ 0) 
and "away" (<^a ~ tt), are typically assumed to contain all 
triggered dijet structure. 

Both the same-side 2D peak and the away-side ID 
peak are broad on azimuth. For minimum-bias jets the 
initial-state parton kt ~ I GeV/c is comparable to the 
minimum-bias mean parton sa 3 GeV/c. The AS 
azimuth r.m.s. width is then approximately 7r/2, and 
the periodic AS peak array is approximated by dipole 
cos (0a) (dashed curve) extending well into the SS re- 
gion [20[. The SS 2D peak for minimum-bias jets in NSD 
p-p collisions is strongly elongated on azimuth with 2:1 
aspect ratio (dash-dotted curve). Thus, the untriggered 
SS and AS jet peaks are strongly overlapping on azimuth. 
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The transverse region invoked in UE studies and indi- 
cated by the hatched regions in the right panel is con- 
ventionally assumed to contain no contribution from a 
triggered high-pt (di)jet (if each jet is confined to a cone 
of radius R < 1) and should therefore be particularly 
sensitive to the UE complementary to the dijet. Fig- 
ure [7] indicates that the TR must include a substantial 
part of the fragment yield from minimum-bias jets. The 
TR is not qualitatively distinguished from any other part 
of the minimum-bias jet structure and may be domi- 
nated by the hard component (minimum-bias jet frag- 
ments) . Since a hard event includes by definition at least 
one jet with mean fragment multiplicity fichj > 2.5 the 
TR in hard events should include a jet fragment density 
(Pjih/ drjdcj) corresponding to at least that multiplicity. 



B. Minimum-bias jets at 200 GeV and 1.8 TeV 

Figure |8] (left panel) shows the surface in Fig. [7] (left 
panel) scaled up by factor 1/0.03 [nj{nch = 2.5) ~ 0.03 
is the mean number of dijets in NSD p-p events within 
A77 = 2] . From spectrum analysis and comparisons with 
pQCD calculations the most-probable jet energy is 3 GeV 
(Q = 6 GeV) As noted, the triangular pair T]- 

acceptance correction has also been removed, restoring 
the accepted jet-related two-particle angular correlations 
in hard p-p collisions, those with a triggered dijet. 




FIG. 8: (Color online) Left: The histogram in Fig. [7] (left 
panel) is renormalized by the calculated jet frequency within 
A77 = 2 from 200 GeV NSD p-p collisions. The pair 
acceptance correction has also been reversed. Right: Projec- 
tion of the histogram in the left panel onto azimuth (solid and 
broken curves) in the format of Fig. 17 of Ref. [l3| compared 
to results from the CDF study (bold horizontal line). 

Figure[5] (right panel) has a format similar to Fig. 17 of 
Ref. 10]. The curves for 200 GeV are for minimum-bias 
jets (within Ar; = 1) with 1.3 (SS) and 1.2 (AS) mean 
fragment multiplicities summing to 2.5, consistent with 
the assumed rich,] in Ref. [2] . The lower solid line assumes 
N^^h = 2.25/3. The soft background is 2x2.5x0.25/50 = 
0.025. Those numbers can be compared with triggered 5 
GeV jets at 1.8 TeV where the computed soft background 
is 2 X 5 X 0.25/50 = 0.05 and the UE hard component 
is (6.9 - 10 X 0.25)/50 « 0.09. The hard component 



inferred from the N±(j)t,trig) plateau value 2.3 is thus 
n/i = 4.4, whereas that inferred from the N± spectrum 
in Fig. [5] is n/i =4 (assuming the same soft component 
for both cases). The two numbers are consistent within 
systematic uncertainties. The 200 GeV estimate only 
establishes a lower limit. Loss of jet fragments at the rj 
acceptance edges for Arj — 1 was not considered in the 
200 GeV spectrum analysis. Such losses may be as much 
as 50% for Arj = 1 and 25% for Arj = 2. A larger fraction 
of partner jets is associated with triggered dijets in the 
larger 77 acceptance. The 5 GeV/c jet trigger at 1.8 TeV 
may also bias jets to larger N± multiplicities compared 
to the untriggered 200 GeV study. Effects of jet bias and 
ri acceptance are discussed further in Sec. IVIIII 



VIII. TRANSVERSE MULTIPLICITY TRENDS 

We can now predict the N±{pt^trig) or <P N ± / drjCLcf) vs 
yt,trig trends observed in UE studies using results from 
200 GeV event-type characterization and angular corre- 
lations. We combine particle density distributions for 
soft and hard events obtained from Eqs. ([5]) and ^ with 
soft and hard spectrum components to determine soft 
and hard event fractions A^; as functions of yt,trig as in 
Eq. ([7]) . To facilitate comparisons with previous analysis 
we adopt the acceptance convention |ry| < 1, pt > 0.5 
GeV/c. Acceptance fractions e for soft and hard spec- 
trum components with a given pt^min lower acceptance 
cutoff and fractions g with given pt^trig upper cutoff are 
determined below. The hard-component contribution to 
the TR has been estimated in the previous section, based 
on angular correlation measurements at 200 GeV. pt,trig 
(single-particle pt, not jet sum Pj) or yt,trig is the trigger 
condition, and the total event multiplicity rich in rj accep- 
tance Ary is the sum of soft rig (projectile dissociation) 
and hard (jet fragment) multiplicities. 



A. pt-cut acceptances 

Figure [9] (left panel) shows the fractions eg and Ch of 
soft and hard particle types that survive the pt accep- 
tance lower bound pt,min or yt,min- The curves are run- 
ning integrals from below of functions So{yt) and 7?o(2/t) 
defined in App. [X] The vertical dashed lines denote the 
CDF acceptance boundary at 0.5 GeV/c that accepts 
about 25% of the soft component and more than 95% of 
the hard component, and the acceptance boundary at 0.2 
GeV/c from Ref. Q that accepted 70% of the pt spec- 
trum (assuming 100% tracking efficiency). 

Figure[9] (right panel) shows factors gs and gh that rep- 
resent the fractions of soft and hard particles surviving 
the single-particle trigger condition pt,trig as an effective 
spectrum upper bound. 
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begins near pt.trig = 2.5 GeV/c, whereas for a leading- 
track-jet trigger the plateau begins near Pj = 5 GeV/c. 



1 1.5 2 



FIG. 9: Left: Running integrals ftom below of soft So and 
hard Ho components defined in App.|X]vspt acceptance lower 
cutoff pt,m,in- The Vertical dashed line indicates a conven- 
tional cutoff at pt = 0.5 GeV/c {yt « 2) for UE studies. 25% 
of the soft component passes that cut. Right: Running in- 
tegrals from above of soft So and hard Ho components vs 
single-particle trigger condition yt,trig that serves as an effec- 
tive upper cutoff on particle momentum. 



B. Transverse multiplicity systematics at 200 GeV 

Transverse multiplicity density cPN^Idrjdcj) is the 
mean 2D angular density within the TR. An expression 
corresponding to Eq. (15) of Ref. [l| is 

iUt.trig) = K{yt,trig)gsiyt.trig)<^snch,s/'2Tr^'n{^) 



dr]d(j) 



K{yt,trtg)gh{yt,trig)e-hn'^hJi/'^T^^'n 



For soft NSD) events at 200 GeV and the CDF yt ac- 
ceptance nch,s/^V ~ 2.5 and « 0.25. For hard events 
n'^f^ J/^T] = 2.25/2 « 1.1 0) and en « 0.95. For sim- 
plicity we assume that the most-probable jets have the 
same properties for any trigger, are not biased. With 
increasing pt,trig the triggered jet energy and fragment 
multiplicity may increase, but the additional higher-pf 
particles should appear closer to the dijet axis and there- 
fore do not contribute significantly to N±_ within the TR. 

Figure [10] (left panel) shows curves from Eq. ([8]) . To 
illustrate the acceptance effect of the trigger condition 
the dotted curve omits both g factors, the dashed curve 
includes gs and the solid curve includes both factors. The 
lower hatched band indicates that only a fraction Cs of 
the soft component passes the pt acceptance. The cor- 
responding density for the full soft component would be 
2.5/27r K, 0.4, larger than the total N± density « 0.3 
predicted for triggered hard events. The upper hatched 
band is dominated by the minijet (most-probable jet) 
contribution to hard events. While the most-probable 
event multiplicity for hard events is 5.7 the soft compo- 
nent is 5.7 — 2.5 = 3.2 just as for soft events. Thus, a 
Pt trigger that selects hard events does not necessarily 
change the soft- component multiplicity. Corresponding 
0.9 TeV CMS data are presented in Fig. 3 (upper left) 



of Ref. [12| . There is reasonable agreement between the 
forms but a significant difference in the plateau values. 
For a leading-track (single-particle) trigger the plateau 
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FIG. 10: Left: Transverse multiplicity angular density vs trig- 
ger condition yt,trig (solid curve) for 200 GeV p-p collisions. 
The single-particle trigger condition also serves as an upper 
acceptance cut on soft and hard components. The dotted 
curve is without that acceptance cut, the dashed curve is with 
the upper cut imposed on the soft component only. Right: 
The same curves plotted on pt,trig- Three labeled hatched 
regions summarize results from higher collision energies for 
comparison. The fourth hatched region estimates an uncer- 
tainty interval for the 200 GeV prediction. 



C. Comparison with other collision energies 

Figure [TU] (right panel) compares the 200 GeV predic- 
tion (solid curve, lower limit) with plateau values from 
higher collision energies. We consider the CDF result 
at 1.8 TeV in more detail based on 2D angular densi- 
ties. The same acceptance pt > 0.5 GeV/c and |r/| < 1 
is applied. For 200 GeV data we have soft compo- 
nent d'^N±,s/dr]d(p = 5 x 0.25/2 x 27r = 0.1 or N±^s = 
5 X 0.25/3 0.4 and hard component d'^ N±,h/dr]d(j) = 
2.25/2 X 27r « 0.18 or N±j, = 2.25/3 = 0.75. The lat- 
ter assumes 2.25/3 jet fragments averaged within the TR. 
For 1.8 TeV data we have soft component d'^N^^s/drjdcf) — 
10 X 0.25/2 X 27r = 0.2 or iV_L,, = 10 x 0.25/3 w 0.8 and 
hard component d'^N±^h/dr]d(j> = 4.4/2 x 27r w 0.35 or 
N±,h4:.A/3 = 1.5. The estimated relatation between 1.8 
TeV and 200 GeV within the TR is about 2:1. 

The calculated 200 GeV TR jet-fragment yield N±^h 
appears to be roughly 50% of the 1.8 TeV yield. However, 
a simulation indicates that for the Arj = 1 acceptance in 
Ref. the mean fragment inefficiency (acceptance-edge 
loss) for minimum-bias jets is about 50%, whereas the in- 
efficiency for A77 = 2 is about 25%. Addition of particles 
outside the rj acceptance but inside the trigger-jet cone as 
in the CDF analysis should not affect dijet (hard) contri- 
butions to the TR. The imposed particle or jet pt trigger 
condition may substantially bias (increase) the jet frag- 
ment multiplicity compared to the untriggered 200 GeV 
study and may bias the soft component to slightly higher 
density (e.g., dus/drj = 5 vs 4). Thus, the 200 GeV soHd 
curves in Fig. [TO] represent a lower limit. The corrected 
plateau value may be substantially higher (illustrated by 
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the hatched region above the sohd curve), possibly com- 
parable to the 0.9 TeV result. 



IX. GLAUBER MODEL, pQCD AND JETS 



that mean participant pathlength v = 2Ni,in/Npart oc 
n^f^ HI]- Those trends in Au-Au can be compared with 
particle- and jet-production systematics in p-p collisions 
where fluctuations may dominate mean values. 



From Sec. I Villi we conclude that imposition of a pt,trig 
condition may select for hard events with at least one jet 
but does not control the soft component, and therefore 
(according to Ref. jH) should not significantly affect the 
p-p collision centrality. Most triggered hard events corre- 
spond to small rich and possibly large impact parameter. 
In Ref. it was observe that imposition of an event- 
multiplicity rich condition controls the jet production rate 
over a large range, and therefore may control p-p colli- 
sion centrality to some extent. We next explore possible 
relations among jet production, imposed rich condition 
and p-p centrality in terms of the Glauber model of A-A 
geometry and pQCD theory. 



A. The Glauber model in A-A collisions 

The Glauber model applied to A-A collisions is based 
on a model nucleon density p^(r), a 3D nuclear density 
function integrating to A nucleons [2l] | . The nuclear pro- 
file function is Ta{s} = / dz pa{z, s), the nucleon areal 
density on a plane spanned by s and perpendicular to 
chord z through the nucleus. The nuclear overlap func- 
tion for intersecting nuclei A and B is 



TAsib) = / d^sTA{s)TBib~s), 



(9) 



a nucleon-pair density on the plane spanned by A-B im- 
pact vector b, with J d^bTAsib) = AB. The number of 
N-N binary collisions with cross section aNN is 



Nbinib) = aNNTAsib), 



(10) 



consistent with basic assumptions of the Glauber model: 
eikonal approximation, straight-line trajectories, lin- 
ear superposition of N-N collisions. The number of 
participant-nucleon pairs in an A-A collision is 

Npart{b)/2 = J d\sTA{s)[l - Po{b ^ s)] < A, (11) 

where Po{b— s) w exp[—aNNTA{b— s)] is the probability 
that a projectile nucleon does not interact in the partner 
nucleus. 

In 200 GeV Au-Au collisions hadron production is ob- 
served to vary with centrality approximately as [l7| 

rich ~ ripp 

{0.9Npart/2 + 0.lNMn) Central Au-Au (12) 
^ npp{0.96Npart/2 + 0MNbin) peripheral Au-Au 

for A?7 — 2. In more- peripheral A-A collisions with N-N 
linear superposition [7| nch{b) oc Npart{b) and jet pro- 
duction scales as rij oc Nbinib) oc ri^l^{b), also implying 



B. A Glauber model for p-p collisions and pQCD 

A similar approach can be applied (formally at least) 
to p-p collisions such that projectile nucleons are treated 
as transverse distributions of partons. Parton distribu- 
tion functions (PDFs) such as fg{x\Q'^) (for gluons) are 
expressed in terms of parton momentum space relative 
to the hadron momentum. To model possible effects of 
p-p centrality we require a differential parton distribu- 
tion function that represents the transverse geometry of 
hadrons. Gluon distribution function fgiXjtlQ"^) is in- 
ferred from DIS measurements {t is momentum transfer 
to the target proton, x is the parton momentum fraction 
within the proton, Q is the parton energy scale). The 
Fourier transform fg{x,s\Q'^) is the inferred gluon trans- 
verse density within a proton, with J d^sfg{x, s\Q^) = 
fg{x\Q^). Gluon density fg{x,s\Q^) is analogous to nu- 
cleon density Ta{s). 

A p-p overlap function can be defined for given 

P2{b)fg{Xl)fg{X2)^ J dhfg{Xl,S)fg{X2,b^m^) 

with J d'^bP2{b) « 1 analogous to Tab (6) defined in 
Eq. The definition assumes an eikonal model with 
straight-line parton trajectories and approximate factor- 
ization of X and s trends within a limited x interval corre- 
sponding to low-energy jets (gluons) near mid-rapidity. It 
also assumes nonfluctuating fs{x,s\Q'^). Functions fg{s) 
and unit-normal P2{b) are assumed to be represented by 
Gaussiaus with widths derived from DIS data 1]. 

The pQCD dijet spectrum for 2-2 parton (gluon) scat- 
tering a, 6 — >■ c, d (collinear factorization) is 



dQ^ 



oc 



J2 [ dxidx2faixi,Q^)fb{x2,Q^){U) 



X daab^cd{xi,X2,Q'^)/dQ'^ 



averaged over the distribution of p-p impact parameters. 
The integral over the dijet spectrum at 200 GeV (with 
spectrum lower bound at 3 GeV) is cTdijet ~ 2.5 mb [gJ. 
We then generalize to s-dependent PDFs fg{x, s\Q'^) and 
obtain for a p-p centrality class with impact parameter b 

d^<Jcd{b,Q^) 



d^bdQ^ 



J dxidx2P2{xuX2,b\Q^) (15) 
a,b 



X daab^cd{xi,X2,Q'^)/dQ'^ 
dacdiQ^) 



P2{b)' 



dQ^ 



analogous to aNNTAB{b) — Nbin{b)- The integral over 
d^b recovers Ec^. p^ . As noted, the second line assumes 
factorization of PDFs within a limited x interval. 
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In the Glauber context rich (dominated by the soft 
component from projectile-nucleon dissociation) may 
serve as a measure of the number of participant partons 
Npart, and observed nj (dijet number) is a measure of 
the number of parton binary colhsions Nhin represented 
by Eqs. ([T4|) and ([T5l) . The question remains whether a 
Glauber description of jet production in p-p colhsions is 
relevant to observations. 



C. Measured jet production vs rich in p-p collisions 

The observed number of minimum-bias (untriggered) 
dijets per 200 GeV NSD p-p collision within acceptance 
Ar/ = 2 is @, [13 



(Jdijet A?7 
r^NSD ^Vi-rr 



0.03, 



(16) 



where A774^ (sa 5 for 200 GeV) is the effective 47r interval 
for a dijet distribution uniform on ry, Udijet ~ 2.5 mb 
is the integral of Eq. (fT4|) [6] and a^sD ~ 36 mb at 200 
GeV [i^l • The NSD result corresponds to an average over 
p-p (or centrality) . If an rich condition is imposed on 
p-p events a p-p centrality dependence may be relevant. 

Fig. [11] (left panel) shows the spectrum hard- 
component multiplicity rih (measured within A77 = 1) in 
the form of ratio rih /n^ vs rich ■ The ratio was estimated 
by two methods, as noted in the panel and described in 
Ref. The observed trend is rih/ris w 0.005 rich/A?/ 
(sohd hue) 0. That result is consistent within ex- 
perimental uncertainties with Eq. ([16]) for NSD colli- 
sions assuming that the mean jet fragment multiplicity 
is TT-chj ~ 2.5 but with A77 = 2. The corresponding dijet 
production vs rich condition is then given by 



rijirich) 



0.03 



{nch/nch.NSoY 
1 + {nch/nch,oY 



(17) 



where nch,NSD/ — 2.5 is assumed at 200 GeV, and 
richfi S> rich,NSD represents the possibility that the 
increasing trend may saturate at some point (dashed 
curve). The number oi jets within Arj is larger by a fac- 
tor representing the fraction of dijets with recoil partner 
in the 77 acceptance, which for Ari/Ari4Tr — 0.4 is about 
1.25 [l3]- Assuming events are Poisson distributed ac- 
cording to dijet multiplicity the probability of soft events 
(no dijet in A77 = 2) is Po{nch) = exp(— n^), the prob- 
ability of hard events (at least one dijet in Arj = 2) is 
1 — Poirich) and the probability of a second dijet (MPI) 
is P2{nch) = [n^jinch) /'^]Po{'nch) (MPI or multiple parton 
interactions is discussed further in Sec. IX Ep . 

Fig. [TT] (right panel) shows the {pt) vs event multi- 
plicity condition for three spectrum analysis methods 
described in Ref. [2]. The "power law" pt spectrum 
model does not describe the spectra well and introduces 
significant systematic errors in inferred {ptj- The two- 
component spectrum model effectively extrapolates the 
spectrum down to pt = 0. In the TCM the fixed mean 
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FIG. 11: Left: Hard-component/soft-component multiplic- 
ity ratio nh/us vs Uch event class for two analysis methods 
(points). The solid curve represents the two-component trend 
nh/us ~ 0.005 Uch/Ar] equivalent to that reported in Ref. [2]. 
The dashed curve estimates possible saturation of the data 
trend. Right: Ensemble-mean pt vs rich event class for three 
analysis methods (points). The curves (described in the text) 
are equivalent to those in the left panel. 



Pt for spectrum soft component S'o(yt) is {pt)s = 0.38 
GeV/c (dashed line). The fixed mean pt for hard com- 
ponent Ho{yt) is {pt)h ~ 1-2 GeV/c. The trend linear 
on Uch (solid line) is then {ns{pt)s + nh{pt)h)/nch as de- 
scribed in Ref. [2], with — Uch — nh and Uh/ris as 
in the left panel. The Uh/ris and (pt) trends are consis- 
tent with minimum-bias jets (minijets) as the underlying 
mechanism for the Uch dependence. Jet production as 
inferred from spectrum and correlation data appears to 
scale as nj oc n^^ « over an extended Uch interval. 
The dashed curve indicates possible saturation of the (pt) 
trend for larger nch, consistent with the left panel. 

Equation (|16p includes the integral Odijct of the par- 
ton spectrum from Eq. ([Ti]) . Equation (|T7)) based on an 
Uch condition may correspond to dadijet/d^b (for some 
6 interval) from Eq. (|15p assuming Uch is strongly cor- 
related with impact parameter h. If such a correlation 
exists parton participant number Npart ~ nch represents 
Eq. pT]) applied to fg{x,s\Q'^). The corresponding num- 
ber of dijets predicted by Eq. (integrated over Q^) 

should be Uj oc Nun ~ n!^Jh according to the eikonal 
assumption adopted by analogy with A-A collisions, for 
which each participant can only interact with target par- 
ticipants along a straight-line trajectory. 



D. Is dijet production related to p-p centrality? 

Is a Glauber model of p-p collisions characterized by 
impact parameter 6, as in Eq. (|15p . consistent with ob- 
servations? We compare s-dependent PDFs from a p-p 
Glauber model in Sec. IIXBI to the observed Uch depen- 
dence of untriggered jet production in Sec. IIX CI We 
assume that soft-component multiplicity Us (~ Uch for 
lower p-p multiplicities) estimates the number of partici- 
pant partons [local parton-hadron duality (LPHD) [23]]. 
The number of participants may he correlated with p-p 
centrality. The number of parton-parton binary encoun- 
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ters is estimated by (proportional to) the observed nuni- 
ber of dijets nj. According to Ref. pj] hard processes 
should arise in more-central p-p collisions. Jet produc- 
tion (probability per event) should then be correlated 
with p-p centrality. Is such a system of rich, jets, parton 
participants, binary collisions, Glauber model and DIS 
data self consistent? 

In general, two extremes are possible: (a) mean-value 
s-differential PDFs derived from DIS data remain fixed 
for all p-p collisions with a given centrality, but different 
p-p centralities (possibly selected by a jet pt trigger) lead 



to different rij cx Nbin and rich numbers per Eq. (|15|) . 
or (b) PDFs fluctuate over a large range, and larger rich 
selects for larger Npart and A^bm with no relation to p-p 
collision centrality. If we accept theory scenario (a) and 

the eikonal assumption we expect rij oc Ntm oc ri^l^ per 
Eq. . If scenario (b) is correct nj oc n^f^ should result 
per Eq. (IT4|) . If rich is not significantly correlated with b 
the system is effectively averaged over b by fluctuations. 
A participant may interact with any participant from the 



collision partner and the rij 



"-ch 



trend follows. 



Figure [T^] represents eikonal scenario (a) (with mean- 
value PDFs, no fluctuations) represented by 



nj{b) 



adijetjb) Ari 
(Tpp{b) A?74^ 

P2(6) 



0.03 X 1.3 



P™(6)' 



where a„ 



In the second line the 
consistent with 



ratio is essentially Nbin{h) ~ ^^part 
the eikonal assumption, and Npart corresponds to fixed 
mean- value PDFs. Figure [K] (left panel) shows functions 
Pa (6) and P„^{h) from Ref. [if (at 7 TeV). Figure [H] (right 
panel) shows the ratio factor appearing in Eq. ([T5|) . The 
factor 1.3 brings the ratio factor to unity for 6 ~ 1.15 fm 
corresponding to NSD p-p collisions (right band). Ac- 
cording to assumptions a jet trigger should select colli- 
sions corresponding to 6 w 0.65 fm (left band) with an 
increase in dijet production by a factor less than 3. 
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Fluctuation-dominated scenario (b) is represented by 

(Tdi]et{nch) A77 



rijirich) 



(rich) Arii 
0.03^^ 
0.03 (rich/^.b)^. 



(19) 



We observe that P2{nch) 



^part 



ilh^iriirich)- In that case 
consistent with an average over 



oc Npart implies that a parton participant 



Nbin 

b, since Nb 

may interact with any target participant. The observed 
nj{nch) trend appears to falsify the eikonal scenario. 

Absolute magnitudes tend to corroborate that con- 
clusion. In An- An collisions Nbin{b) oc N^l^^ propor- 
tional to dijet production varies with centrality approx- 
imately over the interval [1,1000] [l5|. Figure [TT] (right 
panel) shows that for eikonal scenario (a) Nbin{b) varies 
over an interval proportional to [1,3] and cannot describe 
the observed dijet production rijirich) proportional to a 
range greater than [1,100] for an achievable range [1,10] 
of rich/ at 200 GeV. The analogy with A- A collisions 
can be expressed symbolically by 



rij oc Nbin{Npart)P{Npart) 



(20) 



^j^g') In Au-Au collisions with fixed b (or rich) ^ch (x Npartib), 
and P{Npart) ^ 0(1) for a small range of 



N 



bin 



N 



4/3 



part 



Npartib) values (fluctuations are negligible), 
lisions rich ~ Nr, 



In p-p col- 



Vpart but P{Npart) is broadly distributed 
and Nbin oc Npart (fluctuations dominate). 

According to Ref. ,1] by trigging on jets we should se- 
lect more-central p-p collisions with substantially greater 
multiphcities: Nh ^ Ng in Eq. We do observe that 
multiplicity and jet production are highly correlated and 
vary over a large range. If the proton PDF does not fluc- 
tuate where does the large rich in some p-p events come 
from? If rich is strongly correlated with b then large rich 
should correspond to very central p-p collisions. If se- 
lecting dijet events biases to central collisions then a jet 
trigger should result in greatly increased Uch- But the rich 
increase produced by a jet trigger is only a small fraction 
of the available rich range and should correspond to pe- 
ripheral p-p collisions according to the eikonal scenario. 

From the simulation exercise in Sec. I VI I we conclude 
that it is much easier to bias to hard events with a small 
rich increase than with a substantial decrease in 6. From 
the presently-available data we cannot determine the de- 
gree of correlation between rich and b. We cannot rule out 
some degree of correlation for the lowest multiplicities. 
For larger multiplicities correlation seems to be negligi- 
ble and fluctuations dominate. Extension to analysis of 
correlation data vs rich may better access that relation. 



FIG. 12: Left: Probability distributions on p-p impact param- 
eter h from Ref. [3| for dijet production P2 and p-p inelastic 
scattering Pin. Right: Ratio P2/Pin conjectured to be pro- 
portional to the number of dijets per inelastic collision. 



X. DISCUSSION 

This study addresses the relations among pt trigger 
conditions, event multiplicity rich conditions, dijet pro- 



14 



duction and p-p collision centrality. What can be inferred 
from existing data about correlations between pt.trig a-nd 
b, between rich and b, and what is the composition of 
the underlying event? We reconsider the physical mech- 
anisms responsible for observed N±_ vs pt^trig trends and 
N±{pt) spectra. We introduce the concept of a univer- 
sal dijet base that is broad on azimuth and common to 
all jet angular correlations. We conclude that multiple 
parton interactions are unlikely for imposed pt,trig con- 
ditions but may occur with large event multiplicities. 

A. Comparing Nx spectra at different energies 

The collision-energy dependence of 7V_l systematics, in- 
cluding the N±{pt) spectrum, is of considerable interest 
for UE studies. To establish such a trend systematic 
uncertainties and correspondence of analysis conditions 
must be established. A simulation carried out in Secs. lVIl 
and I Villi based on 200 GeV p-p spectrum and correlation 
data is compared to measurements from 0.9 and 1.8 TeV. 
The simulation provides insight into UE comparisons. 

In Fig. [1 the Ni_ spectrum for 0.9 TeV from Fig. 4 
(bottom) of Ref. fl^l (open points) is compared with the 
spectrum for 1.8 TeV from Fig. 37 of Ref. [H (sohd 
points). The rj acceptances differs — Ar] = 4 at 0.9 TeV 
vs A77 = 2 at 1.8 TeV — and the jet trigger condition 
is substantially different — Pj > 5 GeV/c at 1.8 TeV vs 
Pj > 3 GeV/c at 0.9 TeV. Such differences can compli- 
cate interpretation of apparent energy dependence. The 
two-component model can aid such comparisons. 

The soft-component contribution N±,s{pt) at 0.9 TeV 
should be a factor 2 x 3.5/4 « 1.8 times larger than 
that at 1.8 TeV, reflecting differences in 77 acceptance 
and energy-dependent hadron 77 density. The ratio 1.8 is 
consistent with the lowest data points (at 0.75 GeV/c) 
in Fig. [5] Since the accepted soft component is approxi- 
mately one third of A^_l that is a significant issue. 

The hard component N±^h{pt) is influenced by both 
77 acceptance and jet trigger condition. The two-times 
larger rj acceptance at 0.9 TeV should lead to increased 
jet acceptance, but the different trigger condition goes 
in the opposite direction. The higher trigger at 1.8 TeV 
achieves the plateau value 2.3 in Fig. 30 of Ref. [lo| . 
whereas the lower trigger at 0.9 TeV achieves only 50% 
of the plateau value in Fig. 3 (right) of Ref. [l2|- The 
combination may produce comparable hard-component 
yields but a softer hard spectrum shape at 0.9 TeV. 



B. High-pt trigger condition vs p-p centrality 

The trends in Fig.fTUland similar results from UE stud- 
ies such as Ref. can be misinterpreted. According to 
assumptions N± is a "soft" multiplicity not associated 
with the triggered dijet, its variation appearing to match 
the scenario described in Ref. [l]. However, A''^ actually 
includes both soft (nonjet) and hard (trigger-jet-related) 



components. In Sec. I VI B] we demonstrate that while im- 
position of a high-pt trigger condition does bias toward 
hard events (as described by the two-component model) 
it does not significantly increase the soft component of 
N±_ in accepted p-p events, which remains systematically 
consistent with the NSD value. The fraction accepted by 
the effective upper pt.trig cutoff does increase from zero 
at pt,min up a soft plateau value at about 1 GeV/c. 

Given the conjecture in Ref. [l[ the inferred soft- 
component trend implies that the pt trigger does not sig- 
nificantly alter the collision centrality. The event-wise }et 
production probability in peripheral events may be small 
(few percent) but the corresponding fraction of all events 
is large enough (Fig.|3l) that small A^^ (and possibly large 
b) still dominates the overall jet production process. 

The hard-event fraction does increase with pt.trig mitil 
the jet frequency per event rises to unity near pt.trig = 
2.5 GeV/c as described by the TCM. Further increase 
of Pt.trig should bias to higher dijet energy, but the jet- 
related structure (dijet base) within the TR should be 
insensitive to subsequent small-angle jet growth. The 
TCM as applied in the present study with its simplified 
model of the hard component is not intended to describe 
N±{pt) and N±{pt.trig) trends above about 6 GeV/c. 



C. Event multiplicity condition vs p-p centrality 

The relation between Uch and p-p centrality remains 
unclear. In Au-Au collisions Uch and b are strongly corre- 
lated. Fluctuations of geometry and particle production 
play a minor role. In Ref. Q it was established that im- 
posing a condition on Uch in p-p collisions has negligible 
influence on jet properties but directly controls the event- 
wise frequency of minimum-bias dijets. However, in p-p 
collisions fluctuations may dominate mean values, elimi- 
nating any signiflcant correlation between Uch and b ex- 
cept perhaps for peripheral collisions. The low- a; parton 
structure of the nucleon relevant to both minimum-bias 
jets and projectile dissociation may fluctuate strongly. 
Both participant-parton number (cx Uch) and parton bi- 
nary collisions (cx nj) would then fluctuate with some 
broad frequency distribution unrelated to centrality. 

We observe that the two quantities are strongly cor- 
related {uj (X n^fj- Given the form of Eq. that 
relation could arise because the eikonal approximation is 
not valid and all participants can interact or because b 
is effectively averaged over a large interval and the sys- 
tem reverts to Eq. albeit with a large increase in 
the participant number due to fluctuations of the proton 
structure. In the latter case the form of the minimum- 
bias distribution for larger Uch would not relate to the 
p-p impact parameter, only to fluctuations in the proton 
low- a; structure. We conclude that b and Uch might be 
signiflcantly correlated at smaller Uch but may be uncor- 
related at larger Uch- The average b may then retain a 
relatively large value nearly independent of Uch- 
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D. Composition of the Underlying Event 

The event components commonly assumed in analysis 
of the underlying event are (a) an energetic dijet, (b) a 
beam-beam (soft) component representing projectile dis- 
sociation, (c) initial- and final-state QCD radiation, (d) 
possible multiple parton interactions (MPI, additional 
minimum-bias hard scatters) [Tol - [l^ . The UE is identi- 
fied as the combination (b) -|- (c) -I- (d). The transverse 
region on azimuth with multiplicity N± is said to be par- 
ticularly sensitive to the UE: the hard dijet (a) should 
not contribute significantly to the TR. Extrapolation of 
the N± pt spectrum to zero seems to indicate a UE parti- 
cle yield 2.5 x larger than the soft component, suggesting 
a novel production mechanism such as MPI. 

Some of those assumptions can be questioned given 
measured (a) jet fragmentation functions and (b) jet an- 
gular correlations: (a) As shown in Fig. 3 (left panel) of 
Ref. [16l] . with increasing jet energy the fragment den- 
sity (fragmentation function) on rapidity y = ln[(£' + 
p)/mTr] increases at larger y in a self-similar manner. 
Higher-energy jets are built on lower-energy jets by 
adding more higher-momentum fragments which appear 
at smaller angles relative to the dijet (thrust) axis, (b) 
Structure along the dijet axis (y^) can be described as 
a central "string-like" region terminated by "limiting- 
fragmentation" regions at the ends. The central region 
on ijz , its particle density slowly- varying with jet energy, 
corresponds to large relative angles that overlap the TR. 

In Sec. I VIII measured angular correlations for 
minimum-bias jets reveal that no azimuth region excludes 
jet structure. While the 2D distribution is highly struc- 
tured and the SS and AS jet components can be clearly 
distinguished, the projection of minimum-bias jet corre- 
lations onto ID azimuth is almost uniform over 27r, as 
demonstrated in Fig. |8l The jet contribution to the TR 
is then substantial. As the triggered jet momentum or Pj 
increases additional particles at larger pt should appear 
at smaller angles, but the trigger-jet contribution to the 
TR remains the same as for minimum-bias jets. 

Given the results in Sec. IVIII the region \(j>A\ < 0.5 
and |?7a| > 1 would be a better choice to study the UE. 
It is well outside both the SS and AS "jet cones" and 
should exclude most of the central region (on y^) of any 
dijet. Reconstruction of 2D angular correlations relative 
to the trigger-jet axis in future UE studies might provide 
improved (more differential) access to the nonjet UE. 

E. Relevance of multiple parton interactions 

Recent developments in this topic are reviewed in 
Ref. |T3|. Multiple parton interactions have been ex- 
pected to form a major fraction of the UE [l3|- The hard 
component of N±_ inferred from the present study might 
then represent the MPI contribution^ MPI are generated 
in Monte Carlos such as pythia |24| and herwig [111 
conventionally used to model UE systematics. A critical 



aspect of such MCs is the scattered-parton spectrum cut- 
off energy which is treated as an adjustable parameter. 
Significant MPI rates are achieved by lowering the cut- 
off energy so that the effective 2-2 parton cross section is 
comparable to the N-N cross section for a given p-p event 
trigger (e.g., NSD collisions, pt,trig trigger condition). 

A parton spectrum lower bound can be inferred di- 
rectly from p-p and Au-Au spectrum and correlation 
data [i, [i3|. The data in Fig. □ from Ref. ^ indicate 
that the minimum-bias jet frequency for NSD events is 
consistent with a 2.5 mb parton cross section (spectrum 
integral). And the pQCD description in the right panel 
of that figure requires a parton spectrum cutoff near 3 
GeV consistent with (Tdijet ~ 2.5 mb. The present study 
indicates that imposition of a pt trigger condition does 
not significantly bias the soft multiplicity or the collision 
centrality. The probability of a second, minimum-bias 
or semihard dijet accompanying a triggered hard dijet 
should be only a few percent, as for untriggered events. 

If an rich condition is imposed the situation changes 
dramatically since nj oc rv^^ [2|. From Eq. P7|) rich/ ^rj ~ 
13 implies rij « 1, and the probability of at least one MPI 
is « 35%. From Fig. [T] (left panel) we observe that the 
hard component (fragment distribution) does not change 
appreciably with rich, suggesting that the underlying par- 
ton spectrum and 2-2 cross section remain unchanged. 
From the argument in Sec. IIXDI we conclude that the 
rich condition biases to larger parton participant num- 
ber and therefore to multiple parton binary collisions. A 
Pt,trig condition does not significantly alter the (small) 
probability of MPI in p-p collisions. The observed hard 
component of N± is a manifestation of the triggered dijet. 



XI. SUMMARY 

Triggered dijets produced in high-energy p-p collisions 
are by definition accompanied by a complementary un- 
derlying event or UE. The UE is assumed to represent 
several processes, including projectile-proton dissocia- 
tion, initial- and final-state radiation and possibly other 
mechanisms such as multiple parton interactions or MPI. 
It is further assumed that dijet structure is confined to 
regions on {rj, <p) with i? < 1 at and tt azimuth. The 
UE is then studied in a transverse region or TR centered 
at 7r/2 azimuth. The multiplicity N_\_ integrated within 
the azimuth TR and some Ary acceptance is the primary 
UE observable. The pt spectrum and responses to 
imposed conditions are employed to study the UE. 

UE topics of current interest include verifying the phys- 
ical mechanisms that contribute to the TR, especially the 
extent of MPI contributions, and the role of p-p collision 
centrality. The primary observables are N±{pt) spectra, 
triggered N±{pt^trig) systematics and dijet production 
rates. Measured trends have been compared to Monte 
Carlo predictions with a variety of parameter tunes. A 
key issue for Monte Carlos is the invoked pQCD parton 
spectrum model and consequences for the extent of MPI. 
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In the present study we introduce a two-component 
model of minimum-bias and nc/j-dependent spectrum and 
correlation data from p-p collisions at 200 GeV to investi- 
gate UE composition and yield and spectrum trends 
for several collision energies. We test a theoretical con- 
jecture motivated by DIS data that selection of hard col- 
lisions (dijets) with a.pt,trig condition favors more-central 
p-p collisions and increased "soft" multiplicities. 

A two-component analysis of N±{pt) spectra reveals 
that the triggered soft-component yield is not signifi- 
cantly larger than that for minimum-bias collisions, and 
the soft-component spectrum shape may not change sig- 
nificantly over a range of energies. The common hard- 
component shape is consistent with a minimum-bias frag- 
ment distribution predicted by pQCD, and the yields are 
quantitatively comparable with a 200 GeV analysis. 

We observe that a pt.trig condition can effectively select 
for hard collisions, those including a (triggered) dijet, 
but two-component analysis reveals that triggered hard 
events include "soft" hadron densities dus/dr] close to the 
minimum-bias value for a given collision energy. There 
is no indication that such triggered hard p-p events are 
"central" collisions. 

N±{pt,trig) trends are simply explained in the two- 
component context. N± includes soft and hard contribu- 
tions. As noted the soft contribution is consistent with 
(but may be slightly larger than) the measured minimum- 
bias djich/drj. The hard component is consistent with 
minimum-bias jet-related angular correlations associated 
with the triggered dijet itself. The assumption that the 
TR excludes a triggered dijet contribution is incorrect. 
All dijets include a common base corresponding to the 
central region of the dijet on rapidity along the dijet axis 
(string). Those fragments emerge at large angular sepa- 
ration from the dijet axis and contribute to the TR. 

The observed hard component of N±^ might be at- 
tributed to multiple parton interactions or MPI in more- 
central p-p collisions. Larger minimum-bias dijet pro- 
duction can be generated in Monte Carlos by reducing 
the parton spectrum effective lower bound. However, 
parton spectrum properties are strongly constrained by 
minimum-bias p-p data. The parton spectrum integral 
at 200 GeV is approximately 2.5 mb as inferred from the 
rich dependence of hadron spectra and consistent with a 



parton spectrum cutoff neax 3 GeV, implying a minimum- 
bias dijet probability of a few percent. The incidence of 
MPI for an imposed pt,trig condition is small but can 
dominate collisions if a large-rich condition is imposed. 

In conclusion, the two-component model provides a 
good description of the UE as well as other aspects of 
the p-p final state. The status of the theoretical conjec- 
ture that mid-rapidity jet production may depend on p-p 
collision centrality due to the transverse gluon distribu- 
tion in the proton remains unclear. This exercise based 
on yields and spectra points to the need for correlation 
studies, especially the nonjet azimuth quadrupole, that 
may provide additional information about p-p centrality. 

This work was supported in part by the Office of Sci- 
ence of the U.S. DOE under grant DE-FG03-97ER41020. 



Appendix A: TCM model functions 

The unit-integral functions for the two-component 
model (TCM) of rrit or yt spectra as reported in Refs. [1, 
H, 0] are defined here. For 200 GeV p-p collisions the 
soft-component model (Levy distribution on mt) is 

20 

= [l + {m,^m.)/nTY ^^'^ 

with mt — cosh(j/f ), n = 12.8 and T = 0.145 GeV. 
The Gaussian form of the hard-component model is 

Hoiyt) = 0.33exp{-[(yt - yto)/cTyf/2} (A2) 

on yt, with yto — 2.67 {pt ~ 1 GcV/c) and dy^ = 0.445. 
The coefficients (determined by the unit-integral condi- 
tion) depend on the specific model parameters. In Ref. [1] 
the hard-component model function is generalized to a 
Gaussian with power-law tail to accommodate the under- 
lying parton energy spectrum. In Ref. ^6] the calculated 
pQCD hard component derived from fragmentation func- 
tions that describes the spectrum data accurately also de- 
viates from a Gaussian below the distribution mode. The 
simplified Gaussian form adopted for the present analy- 
sis is compared with the 200 GeV pQCD calculation and 
p-p data in Fig. [1] (right panel). 
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